Denaturation of Circular DNA: Supercoil Mechanism 
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The denaturation transition which takes place in circular DNA is analyzed by extending the 
Poland-Scheraga model to include the winding degrees of freedom. We consider the case of a 
homopolymer whereby the winding number of the double stranded helix, released by a loop de- 
naturation, is absorbed by supercoils. We find that as in the case of linear DNA, the order of the 
transition is determined by the loop exponent c. However the first order transition displayed by 
the PS model for c > 2 in linear DNA is replaced by a continuous transition with arbitrarily high 
order as c approaches 2, while the second-order transition found in the linear case in the regime 
1 < c < 2 disappears. In addition, our analysis reveals that melting under fixed linking number 
is a condensation transition, where the condensate is a macroscopic loop which appears above the 
critical temperature. 

PACS numbers: 87.15.Zg, 36.20.Ey 



I. INTRODUCTION 

Thermal denaturation of DNA is a process by which 
the two strands of the molecule unbind upon heating. A 
good understanding of the underlying physics is relevant 
to certain biological systems (e.g., thermophilic organ- 
isms [J, @|) as well as synthetic technologies Q such as 
polymerase chain reaction (PCR) [3, [j| and DNA mi- 
croarrays @. The unbinding transition takes place at 
a specific temperature, coined melting or denaturation 
temperature, which can be defined experimentally as the 
temperature at which the fraction of unbound base pairs 
reaches, say, half of its maximal value. For a relatively 
homogenous DNA chain composed largely of A-T (or 
G-C) pairs, melting takes place through a very sharp 
increase in the fraction of broken bases, suggesting a 
first-order phase transition in an idealized homogeneous 
system. This phase transition has been investigated by 
means of various theoretical approaches developed in re- 
cent decades 0-Qil]- 

A prototypical model employed in theoretical stud- 
ies of this phenomenon is the Poland-Scheraga (PS) 
model [9j in which a microscopic configuration of the 
DNA molecule is described by an alternating succession 
of bound segments (dsDNA) and denaturated loops (ss- 
DNA). As the temperature is increased the total length 
of the bound segments decreases, eventually vanishing at 
the melting transition. The transition is a result of the 
competition between the enthalpy associated with the hy- 
drogen bonding of the matching bases, and the entropy of 
loops. The loop entropy has the asymptotic form ~ s l /l c 
for large loop size I, where s is a geometric, non-universal 
constant and c is a universal exponent. The original PS 
model makes the simplifying assumption that the bind- 
ing energy is the same for all base pairs, in which case the 
nature of the transition depends only on the parameter 
c [9j. For c < 1 no transition takes place and the two 



strands are bound at all temperatures. For 1 < c < 2 the 
model exhibits a second-order melting transition where 
the average loop length increases and becomes macro- 
scopic of order L as the critical point is approached from 
below. For c > 2 the transition is first order and the av- 
erage loop length remains O(l) for T < T c . For T > T c a 
macroscopic loop, formed abruptly at T c , is present. In 
d = 3 dimensions and with exclusion interaction prop- 
erly taken into account, one obtains c « 2.12 (Til. [l9L |2Q] 
and the transition is predicted to be first order. The PS 
model has later been extended to address the sequence 
dependence of the melting transition in heteropolymeric 
DNAs HU- 

The DNA molecule is helical, and therefore denatu- 
ration entails unwinding of the two strands around one 
another. The PS model ignores this fact, as the elas- 
tic strain can be relaxed by the rotation of the chain 
ends. However there are cases where the helicity can not 
be ignored. For example, bacteria have circular DNAs 
(plasmids) whose linking number (the number of times 
one strand winds around the other) is a topological in- 
variant. Similarly, certain single-molecule experiments 
require the chain ends to be rotationally constrained. In 
such cases, unwinding of a loop is possible only if some 
extra linking number can be absorbed by the rest of the 
molecule. 

Previous studies that model denaturation of circular 
DNA proposed two mechanisms by which bound DNA 
segments may host extra linking number released by 
opening a loop: (a) increasing the twist (the excess stack- 
ing angle integrated along the centerline) [22|, [23| ; or (b) 
increasing the writhe (which is a function of the centerline 
configuration itself), for example by forming a supercoil 
[H-[26j]. The AFM images of thermally denatured DNA 
circles adsorbed on a mica surface suggest that supercoils 
do form in conjunction with denaturation loops [27| . Nu- 
merical studies similarly point at the writhe as the dom- 
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inant mechanism for absorbing the extra linking number 
in long DNA circles [26j . 

In this paper we study in detail the case of supercoils. 
In an earlier work this model has been studied at tem- 
peratures below the melting point [24j . by means of a 
grand canonical treatment where the expectation value 
of the linking number is fixed. Here we generalize this 
approach and further consider the high-temperature de- 
natured phase in order to study the nature of the melting 
transition. The validity of our results is then verified by 
a direct calculation within a canonical formalism where 
the linking number is strictly conserved. This approach 
allows us to point out an inconsistency in the assumed 
analogy with the PS model in Ref.[24|. Finally, we find 
the following phase diagram: For c < 2 the model ex- 
hibits no phase transition and a steady increase of loop 
fraction with temperature. For c > 2 a continuous tran- 



takes place, where \q] is the upper 
The order of the transition tends to 



sition of order 

integer value of q. 
infinity as c — ¥ 2. 

The paper is arranged as follows: In section [Til we 
present the model. In section Hill the denaturation tran- 
sition is first established in the grand-canonical ensemble, 
where we introduce a regularization scheme used earlier 
in [l9| . This procedure allows us to draw an analogy be- 
tween the high-temperature phase and a Bose-Einstein 
condensate where a critical fluid (microscopic loops) co- 
exists with a condensate (a single macroscopic loop). In 
section ITVl we reinvestigate the model within the canon- 
ical formalism: while we observe a general agreement be- 
tween the two ensembles, we also point out a difference 
between the corresponding condensates that suggests the 
inequivalence of the two ensembles for finite systems in 
the present context. Finally, in section [Vj we present 
some concluding remarks and discuss possible future di- 
rections. 
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FIG. 1: A typical configuration of a portion of the circular 
DNA model used in this study. 



"cooperativity parameter" (3. The (universal) 
loop exponent c ~ 2.12 is determined by the di- 
mensionality of the embedding space (= 3) and the 
connective topology of the polymer system [ll], HH ■ 



3. a supercoil, in which two halves of a dsDNA seg- 
ment wind around each other (see FigfT]). The cor- 
responding Boltzmann weight is given by e~ l ^ Es = 
v , where E s (0 > E s > E^) is the energy gain of 
a base pair in a supercoiled segment. Our model 
reduces to the PS model when v = 0. 

It is assumed that supercoils occur within bound regions 
only and hence a loop is always terminated by two bound 
segments (of type 1 above). A typical configuration of 
part of a circular DNA molecule is shown in FigfT] where 
l\ denotes the length of the i th loop, while l\^ and Ifj 
stand for the lengths of the j th bound segment and the 
j th supercoil following the i th loop, respectively. The 
Boltzmann weight corresponding to the configuration in 
Figfflis 



II. MODEL DEFINITION 

Following [iH we extend the PS model to include su- 
percoiled DNA segments. Thus, a microscopic configura- 
tion is composed of an alternating arrangement of three 
types of segments: 

1. a bound segment, in which base pairs are intact 
but no supercoiling takes place. Following the PS 
model, we neglect the entropic contribution of such 
a segment, so that its Boltzmann weight is solely 
determined by the binding energy Et> < and the 
segment length I as e~^ lEb = w l , where (3 = l/fe^T. 

2. a loop, in which pairing is sacrificed in favor of en- 
tropy as the persistence length of ssDNA is roughly 
10 times shorter than that of dsDNA. The asso- 
ciated Boltzmann weight is purely entropic and 
asymptotically given as O (Z ) = Aj^, where s is a 
geometrical factor and A is a constant coined the 



fl x u^-i.o X fi (Z?) X cA° X l/lo X 

aA 1 x v l U x u 1 ^ 2 x n (z' +1 ) 

Let Lb, L s and L/ be the total length of bound, super- 
coil and loop segments respectively. The length of the 
DNA is given by L& + L s + Li = L. The conservation of 
the linking number is imposed by the additional condi- 
tion that an increase in the total loop length Li (reducing 
the linking number) is compensated by a proportional 
increase in the total supercoil length L s (recovering the 
linking number) , and vice versa. Given the ground state 
Lb = L, L s = L; = 0, this yields the constraint L s = aLi. 
a is the proportionality constant and for simplicity we 
assume here a — 1, though the result is qualitatively 
the same for other values [25|]. In this model, as in the 
PS case, it is more convenient to work within a grand 
canonical ensemble, where the above constaint is relaxed 
to the equality of corresponding ensemble averages, i.e., 
(L s ) = {L l ). ' 



3 



III. GRAND CANONICAL TREATMENT 

For completeness we first outline the derivation in [24| 
for this case. To account for the two constraints above, 
the grand partition sum is constructed as a function of 
two fugacities z and [i as 



Q(z,fi) 



Z(Lb, Li - L s )z ijl 



L . .L, —L. 



(1) 



where Z(L b , L s — L\) is the canonical partition sum. Note 
that for fj,= l, Eq. ([1} is the grand canonical partition 
function of the Poland-Scheraga model extended to in- 
clude all possible supercoil segment insertions. While 
this partition sum is different from that of the original 
PS model, it qualitatively yields the same phase dia- 
gram [24j |. 

The values of z and [L are set by the conditions 



L = ^ (=L b + L s+Ll ) 







d log z 
<91og<2 
<91og^ 



(= U - L s 



(2) 
(3) 



Assuming that there is at least one bounded base pair, 
the grand partition sum can be written as 



Q{z,n) = V(z,[i) + V(z,[i)U(z[i)V(z,n) + 
= V(z,[i) 

l-V(z,[i)U(z[xy 



(4) 
(5) 



with 



V(z,[i) 



V(z) 



1 - v(z)w{z/fj,y 

00 

z — ' 1 — Ul. 



w{z/n) = £ U~ 



00 / \ n 

z \ vz 



* \ [I J [I — vz 



(6) 
(7) 

(8) 



(sz/z) T 



A<S> c (sz[i). (9) 



The functions f7, V and represent the grand partition 
sums for loops, bound segments and supercoils respec- 
tively. The polylog function $ c (<z) is given by 



*c(?) = E 



n=l 



(10) 



It is an analytic function everywhere except for a branch 
cut at q € [1, 00) [2^|. It satisfies the relation 



^-* c (g) = -* c -i(g) 

dq q 



(11) 



By inserting Eqs. (JSjlQj) into ([5]), Q(z,[i) can be written 
as 



Q(z,n) 



1 



1 



(12) 



From this explicit form the constraints given by Eqs.([2][3]) 
are readily transformed into 



1 

LOZ 



1 



vz 



fl — vz 
vz 
(p — vz) 2 



= A$> c (sz[i), (13) 
= -$ c _r( S z M ), (14) 



where z and [i from here on refer to the corresponding 
values in the thermodynamic limit (L — > 00) which is 
assumed in the derivation of Eq. p^l) . Denoting by nit, 
m s and m; the average density of base pairs in bound 
segments, supercoils and loops, respectively, one finds 
that 



m b = 



mi 



dlogz 
d log ui ' 
d log z 
d log v ' 
d log z 
d log s ' 



(15) 
(16) 
(17) 



It is more convenient to work with the transformed vari- 
ables x = szfJ, and y — vz/ [i. Physically (x/s) is the 
fugacity associated with a unit increase in the total loop 
length, and (y/v) is the similar fugacity of supercoils. 
Under this change of variables Eqs. ([TBI [T4"|) become 



yfsv 1 



wy^ 1-y 
V 



= v4$ c _i(x) 



(i - y) 2 

Considering y as a function of x through Eq. (fT9"f , let 



(18) 
(19) 



A$ c (x) 



so that Eq. (fl8|) can be written as 



y 



(20) 



G{x) 



j(E b -±E s ) 



H(T). (21) 



Note that y and G(x) are increasing functions of x in 
the physically relevant regime < x < 1. The lower 
bound x — is achieved at zero temperature since 
Ym\T^o H{T) = 0, while the upper bound is unity since 
$ c (x) is a divergent sum for x > 1. The presence 
of a thermodynamic phase transition then depends on 
whether x(T c ) — 1 is achieved for some finite tempera- 
ture T c . Two regimes emerge as shown in Fig[2j 
(i) c < 2: lim x _yi $ c _i(x) = 00, therefore Eqs. (|18ll9p 
have a solution in the interval < x < 1 at all tempera- 
tures, 
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FIG. 2: A plot of G{x), as defined in (f20"|) . as function of 
x. The solid line corresponds to c = 1.5 while the dashed 
line to c = 2.5. The horizontal lines corresponds to different 
temperatures. The arrows point to the solutions of Eg. (|18p 
which is equivalent to the thermodynamic limit. While for 
c = 1.5 there exists a solution for all temperatures, for c = 2.5 
there is a solution only up to some finite temperature, in this 
case v^ 2 uiZ X ~ 0.25. The parameters which were used in this 
plot are E b = -3, E a = -2, s = 5 and A = 0.1. 



(ii) c> 2: Note that $ Q (1) = J2 n 1 / n ° = Ca is the Rie- 
mann Zeta function [3(J. Then, as x — > 1, the RHS of 
Eqs. (|18ll9[) remain finite since oo > f c _i > ( c > 1. For 
suitable values of A and s (which guarantee G(l) < 1), 
there exists a temperature T c such that x(T > T c ) = 1. 
The resulting nonanalyticity at T c translates into the sin- 
gular behavior of other quantities like the density rrib 
which underlies the melting transition. T c is given by 



G(l) 



Below T c , the system is fully defined by Eqs.(fT 
Above T c one has 



SZfJ, 



1 



(22) 



(23) 



and an additional equation is necessary to impose the 
two constraints above. To this end, we follow Ref.[l9j 
and introduce a cutoff M on the maximal loop size. In 
this reduced ensemble the partition function is analytic, 
thus Eqs. (fF?iri4"l) are valid at all temperatures. The limit 
AI — y oo reveals precisely how these equations are modi- 
fied above T c , as discussed below. 



A. Regularizing the Grand Canonical Ensemble 

Introducing an upper cutoff M on the allowed loop 
size, the loop partition sum U (z/i) is replaced by 



71=1 



where <&^r (?) is the "loop-truncated" Polylog function, 
while the relation ■^Qf(sq) = | < I ) e-i( s <?) still holds. 
The grand canonical partition sum is then 



0,m) = 



1 - V(z,fx)U M (z^) ' 
and Eqs. (|2l3p for the constraints can be written as 



— -1 -- 

OJZ J fi- 



vz 



A$f{szpL), 
A_ 



(H - vz) 2 



(24) 
(25) 



These equations hold for all T, since <&^r (x) is an ana- 
lytic function. Our goal now is to analyze these equa- 
tions in the limit M — > oo for temperatures above T c . 
In this approach one should, in fact, consider the grand 
canonical ensemble with a finite but large average length 
of the DNA molecule (L). One should then consider 
the limit M, (L) -> oo with M - 0((L)). While con- 
sidering finite (L), Eg. (|2~41 is no longer exact but has 

({L) X ^ correction. However, this correction does not 

modify the analysis presented below as it vanishes in the 
limit M ~ (L) —> oo. We therefore take (L) = oo and 
then M -> oo. 

Let T C M be the temperature at which szfi = 1 in this 
loop-truncated model, so that for T > T^ 1 we have szfi — 

1 + e(M,T) with e > 0. Clearly, as M -> oo, T* 1 -> T c 
and e — > so that sz[i = 1 for all T > T c . Then, for a 
given temperature T > supm (Tc ' 1 ) we have 



1 

OJZ 



fj, 



fi — vz 
vz 

(/i - vz) 2 



A{$M(l) + b(e)), 
1 (l)+a(e)) 



where 



szfi = 1 + e, 



(26) 
(27) 

(28) 



and where a(e) and b(e) are cutoff-dependent corrections 
at least one of which is nonzero (otherwise the system is 
overdetermined) . We continue by assuming that Me — > 
oo as M — > oo, and checking that this assumption is self 
consistent. With this assumption the leading behavior of 
a(e) and 6(e) is found as 



A I 



(i 



i 



71=1 

i-M ex 
dx- 



rv 



x'- 



(M 



and 



M 



He) - E 



(1 + e)' 



a(e)/M, 



(29) 



(30) 
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Therefore, the only cutoff-independent choice is a(0) = 
ao and 6(0) = for some constant a . Moreover, the 
asymptotic form of e as a function of M follows from 
Eq.d^J) as 

e( M,r) = (c-2>^ + o(™V (31) 



M 



M 



demonstrating the self-consistency of the assumption 
above (see also HH). We conclude that while for T <T C 
Eqs. (|13lll4")) hold, above T c they are replaced by 



1 



U!Z /i — VZ 

vz 



A 



(32) 



, 2 - — (Cc-i +a ), (33) 
(/i - ^) M 



sz/i = 1, 
from which we can extract ao 

v 

a 



As (/i — vz) 



Cc-l 



(34) 



(35) 



It measures the density of base pairs that reside within 
a macroscopic loop - or condensate - that appears above 
T c . In the next section we discuss the order of the phase 
transition, where we take into account the condensate 
correction which was omitted in Ref. I2JI 



Below T c this can be obtained by noting that Eq.(|14p 
implies §^ = 0. Equation (f31)f can then be used to cal- 
culate t^j and and finally Eq. ([T4")l is used again to 
eliminate the polylog function. Above T c Ea.(|14p docs 
not hold and is replaced by szpL = 1. Equation (|39[) is 
then obtained by evaluating the partial derivatives ap- 
pearing in Eq. (|38|) . Equation (|39|) implies that the order 
parameter m& is continuous across the transition since z 
and (j, are continuous functions of the temperature. Thus 
the transition is continuous. 

For a detailed analysis of the singularity it is conve- 
nient to express mb in terms of the x,y variables. Let 
Sx, Sy, and Smb denote the deviation of x, y, and mi,, re- 
spectively, from their values at T c due to a slight change 
in temperature t = T — T c . First we explore the relations 
among t, Sx and Sy. Above the transition x = 1 and 
hence Sx = 0. Thus & c (x) in Eq. (fT8|) becomes £ c - As a 
result y has a power series expansion above T c where to 
leading order Sy oc t. Hence for t > 0, Smb = f(t) where 
/ is analytic near t = 0. Below the transition one has to 
make use of the expansion of the polylog function 



$ c _!(l-fe) = Cc-i+Cc-2Sx+...+T{c-2)5x c 



(41) 



where T(c— 2) is the Gamma function and the last term is 
the leading singular term in the expansion. We proceed 
by separately considering two regimes of the parameter 
c. 



B. Order of the Transition 

We show below that the above phase transition is con- 
tinuous and then investigate the nature of the singularity 
at T c . Consider the fraction of base pairs in bound seg- 
ments, i.e., mj, = —dlog(z) / dlog(u). Defining 



It 



-A^> c {szn), (36) 



and noting that P{z 1 h;uj) = for the poles of the par- 
tition function, the smallest of which yields the thermo- 
dynamic limit, we find 







dP 

duj 



dP 

duj 



dP dz 
dz dui 



dPdfj, 



Rearranging and using Eq. (|15l) 

ui dPj duj 



m b 



z dP/dz + (dP/dfi) (dfi/dz) ' 



(37) 



(38) 



Evaluating the derivatives, making use of Eqs. (|13ll4[ 
I23I36P we find that both below and above the critical 
temperature to& is given by 



• For 2 < c < 3, the expansion (|4Tj) becomes 
$c_i(l - Sx) w Cc-i + T(c - 2)Sx c ~ 2 and there- 
fore Eq-lflg]) yields 



Sy ~ Sx c 



(42) 



which implies Sx <C Sy. Thus in the vicinity 
of the transition temperature Eqs. (fT8l [40|l yield 
Smb ~ f(t) + aSx where f(t) is the same function 
as above the transition, and a is a constant. Using 
(I42p and noting that Sy oc t one finally obtains 

Smb — f(t) + at^ 1 where a is a constant. Since 



f(t) is analytic function, the 



derivative of 



mb is discontinuous and the transition is of order 



• For c > 3, the expansion of the polylog function 
is $ c -i(l - Sx) w Cc-i + (c-2Sx. Hence Eq.lfigp 
yields Sy oc Sx, which together with (1181) implies 
Sx oc Sy oc t. Thus Smb ~ f(t) + yt where 7 is 
a constant. This implies that the first derivative 
of Smb is discontinuous, and the transition is of 
second order. 



m b = 



2ivu/j,z 2 
(fx - vz) 

luyJxy 



sv(l - y) 



(39) 
(40) 



In summary, the transition is characterized by the sin- 
gular behavior of Smb below: 



Smb = 



f(t) 

/(*)+«*" 



t > 
t > 



(43) 
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with 



1 

c-2 
1 



2 < c < 3 
c > 3 



(44) 



where f(t) can be expressed as a power series in t for 
t > 0. Since to/ = m s = (1 — to,{,)/2, a similar singular 
behavior is exhibited by these variables. Hence the de- 
naturation transition of a circular DNA is second order 
for c > 3, third order for 2.5 < c < 3, forth order for 
| < c < 2.5, etc., approaching infinite order as c — > 2. 
No phase transition takes place for c < 2. In contrast, 
a DNA without helicity (as described by the original PS 
model) melts through a first order transition for c > 2 
and a second order transition for 1 < c < 2. 



C. High Temperature Phase 

The high-temperature phase of the PS model is com- 
posed of an all-encompassing macroscopic loop created at 
T c through a jump in the loop fraction to its maximum 
value mi = 1. Here, we not only have a smoother tran- 
sition but also a qualitatively different denatured phase. 
For example, the loop fraction reaches its maximum value 
(to; = 1/2 within the present model) only as T — > oo and 
it continuously increases across and above T c . At this 
point, one is tempted to ask what has changed qualita- 
tively across the transition. In this section, we show that 
a macroscopic loop is again the distinguishing feature. 
However, instead of being an all-or-none phenomenon, 
the dominance of the macro-loop among the denatured 
base pairs grows steadily from T c on. Below we analyze 
the loop length distribution, Pm(1), demonstrating that 
in addition to the power law behavior on microscopic 
scale, it exhibits a peak at lengths of order M whose in- 
tegrated weight is of order 1/M. This peak represents 
the macroscopic loop which opens up above T c . Note 
that the probability distribution functions for bound and 
supercoiled segment lengths are still exponential in the 
length n, since the corresponding Boltzmann weights are 
(ujz) n and y n , respectively. 

The loop size distribution in the "loop-truncated" 
model is given by 



Pm(1) 



1 



(szfi) 1 



<f>M(szfi) l c 



&(M-l), 



(45) 



where O(x) is the Heaviside function. Differentiating 
with respect to I and noting that szfi — e e we find that 
this distribution exhibits a minimum at 



r 



M 



c - 2 log M ' 



(46) 



where Eq.(f3"Tj) has been used. The large I distribution is 
peaked at I = M with p M (M) Q l e eM /M c ~ M~ 2 . 
Thus the integrated weight of the peak is 0(1/ M) up 
to logarithmic corrections. Since as discussed above 
M ~ O(L), One expects 0(1) number of macroscopic 
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FIG. 3: The loop size distribution p(l) in the canonical (solid 
line) and the regularized grand-canonical (dashed line) ensem- 
bles. For small values of I the critical phase can be identified 
(where p(l) ~ l~ c )- For I ~ L the canonical curve shows the 
"bump" around I ~ £ (see text) while the grand-canonical 
curve behaves in a somewhat different manner. The param- 
eters which were used for this plot are c = 3.5, L = 400, 
M = 200, E b = -3, E s = -2, s = 5 and A = 0.1 and T — 3 
{T c = 1.167) 



loops to open up above T c . In fact one can argue that 
entropy favors a single macroscopic loop [3l| . To see this 
one can compare the probability of a state with only one 
macroscopic loop with that of configurations with two 
macroscopic loops: Assuming that there are l con ~ M 
base pairs within the condensed phase, the weight of con- 
figurations with single loop is 



Tl (Icon) ~ LlPM (Icon) ~ O (L 1 c ) (SZ/J,) 



The weight of configurations with two macroscopic loops 
is 



r 2 {ic 



L 



^2 p m ( n ) p m (' c °« _ n ) 



O (L 3 - 2c ) (szn) L -> 



As c > 2, it follows that configurations with a single 
macroscopic loop dominates the ensemble in the limit 
M ~ L -> co. 

The condensation phenomenon observed in this model 
is reminiscent of condensation in Bose-Einstein Gas and 
the zero-range process (ZRP) when the density of par- 
ticles is above a critical value [HI, HI] ■ Figure [3] shows 
the loop size distribution for finite M and T > T c . A 
power-law decay with the exponent c for I -C M and a 
peak for I < M which is the precursor of the 5-function 
representing the macroscopic loop are evident. 



IV. CANONICAL TREATMENT 

In order to justify the regularization procedure ap- 
plied in the grand canonical ensemble we study the model 
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within the canonical ensemble, namely with fixed L and 
L; — L s = 0. In addition this approach allows us to study 
the properties of the condensate and to further illuminate 
the mathematical structure underlying the phase transi- 
tion. The canonical partition function can be obtained 
from the grand sum in Eq. dTJ by means of Cauchy inte- 
gration: 



Z(L,Li-L s ) 



1 



dfi 



dz- 



(2n) z / c w ' Jew z i +V L '- Ls+1 ' 

(47) 

where and C^ z ' are circular, counter-clockwise ori- 
ented contours which are centered at the origin and en- 
close no singularity of Q(z, /i) (Figj4j) . Enforcing the link- 
ing number constraint, Li — L s , and using Eq.(|12p yield 



Z(L,0) = 



I{z,n) 



l 



1 



A<Z> c (szn) 



yL + l 



(48) 



(49) 



For \n\ sufficiently small so that \sz/j,\ < 1, let zo be the 
nontrivial pole of I{z, fi) in the z-plane, given by Eg. (1131) . 
Then, by Cauchy's integral theorem, the integration con- 
tour can be replaced by + C^} shown in FigJU 
Due to the factor z~ L in (|49p the dominant contribution 



comes from G 



0) 



and we obtain 

J- I *>■ 

2m 



Z(L,0)~-^f m^^-i-i. (50) 

We now evaluate the integral separately below and 
above the critical point. Below the transition, the inte- 
grand in Eq. (l50)) has a saddle point given by 7^7 | w = 
and \no\ < l/s\zo\. The partition function can now be 
evaluated by first deforming C'^ into the contour Cg 
which passes through this saddle point (Fig(4j) and then 
approximating the integral by the contribution from the 
vicinity of fj,Q, i.e., 

-Llog z (/io) 



Z(L, 0) ~ e 



(51) 

After differentiating Ea.([T3| with respect to /i and set- 
ting ^ = we find Eq.fTJJ as the saddle-point condi- 
tion. These two equations fix z and /io and describe the 
system for T < T c , as was found earlier in the grand- 
canonical framework. Note that the free energy is ob- 
tained from z through Eq. (j5Tj) . 

Above the critical temperature, this procedure is not 
applicable, as the solution of = for (1q now lies 
on the branch cut. However, it is found that Eq. (|5"Tj) 
holds, with zq and /io given now by Eqs. (|13l23l) rather 
than (|13I14I) as obtained within the grand canonical en- 
semble. This can be shown by evaluating the integral 
in Eq. (f5T))) along another contour C%£ shown in Fig 0] on 
which |szo^| < 1. After a change of variables e~ u — szo/x, 
Eq. (l50l) transforms to 



Z(L,0) 



Llogzo(ti) 



(52) 



where a(u) — [\ + d j is a nonextensive correc- 
tion to the free energy that can be neglected. The main 
contribution along the contour Cj£ is from the neighbor- 
hood of the branch cut where \u\ <C 1 with Re [it] positive 
and as small as desired. We therefore express log zq(u) 
in terms of the small parameter u by using the implicit 
equation (|13p and the nonanalytic expansion of $ c (l — u) 
given by Eq. (l4"T1) . to obtain 



log z (u) 



Lc-ij 

J2 h ^ r ' 



b c u c 



where b n are temperature dependent coefficients with 
b = logz (0) and b n = ±J^p>\ u=0 . The coefficient of 
the linear term b\ vanishes at T c . This follows directly 
from Eqs. (|13ll4l23p . It changes sign from b\ > be- 
low the transition, where szoMo < 1, to 61 < above 
it. Let g(u) be the nonlinear part of the expansion 

~ Eil^ ^ + b c U C 

g(u) — b c u c ~ l . One therefore has 



En= 2 1J b nU n + b c u c - 1 . Note that for c < 3, 



Z{L,0) ~ e -£l°«*o(0) 



du 
2ttI 



-L[biu+g(u)] 



(53) 



As u is approximately imaginary in the region of inter- 
est, the integrand is oscillatory, yielding vanishing con- 
tribution at large L except in the small region where 
lm[u] < 0(7;)- As a result one may expand the in- 



tegrand in Eq. (|53|) as e" 



-L[b 1U +g(u)] 



-bi Lu 



[l-Lg{v)]. 



Moreover, the integration contour can be replaced by the 
right vertical tangent of in Fig|4j Combining these 
observations we get 



Z(L,0) ~ e - Llog2o(0) 



du e 



-biLu 



[l-Lg{u)\. (54) 



The analytic terms of the integrand do not contribute, 
since the integration yields a delta function 5 (b\L) or 
its derivatives [32]. Therefore, the partition function is 
determined solely by the nonanalytic term in g{u) as: 



Z(L,0) 



e -Llog* (0) b t 



-Llog z (0)_ 



due-^ Lu u c - x 



c r.c-i ' 



b\L 



(55) 



where b = b c ^^T(c) The free energy density 

is, of course, continuous across T c and above the critical 
temperature it is determined by Eqs. (|13l) and (|23p . as in 
the grand canonical treatment. 



A. High Temperature Phase 

In this subsection we consider the loop size distribution 
p(l) at temperatures above T c . As in the grand canoni- 
cal ensemble, a condensate phase composed of a macro- 
scopic loop is found, although the details of the peak in 
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(a) (b) 
Im [z] 




FIG. 4: The integration procedure used for the canonical partition function: (a) In the z plane, the contour C*- 2 ' which encircles 
the origin can be replaced by Cp Z ^ around the pole at zo and which wraps the branch cut (thick line) and closes at infinity; 
(b) in the fj, plane, the contour C^- 1 can be deformed to pass through a saddle point [iq when it exits. Otherwise the dominant 
contribution comes from the vicinity of the branch point (see text). 



p(l) corresponding to this phase are different. The anal- 
ysis follows the analysis carried out for the condensation 
transition in the zero-range process [32| . Here we just 
outline the main results. 

Within the canonical ensemble the loop size distribu- 
tion is given by 



p(l) = A 



i l Z{L-l,-l) 



l c Z(L,0) ' 
where Z(L — I, — I) is given by 

I f , f , Q(z,fx) 



(56) 



Z(L - 1,-1) 



(2ni) Jew JC* 



dz 



-i 



2ni 



cm 



with 



du e _ - L [ logz °( u ) + ^ u ] 

27TJ JqM 

l/L and e~ u = szq(i. Expanding for small u 



yields 



P(l) 



AI(l/L) 
T c 7(0) ' 



(57) 



I{4>) = — / dwe -£[(&i+0)«+s(«)], (58) 



For T > T c , the function I(4>) develops a peak at <f> ~ 
— b\ = £. This is demonstrated separately for 2 < c < 3 
and c > 3. 

(a) For 2 < c < 3, the leading-order term in g(u) is 
the nonanalytic term u c and I((f>) can be written in the 
form 



1(0 =Z-V(°-i)y c 



(59) 



The asymptotic behavior of the scaling function V c [q] are 
given by 



Vc[q] 



a\q\ q 

7 (3-c)/2(c- 2) p -c 2 g< t! - 1 V 



-00 



ciq y 



q 00, 



(60) 
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where the constants a, c\ and C2 are given in Eqs.(81- 
83) of 32]. EquationsjlTlinni) together with Eq.([55) yield 
after some algebra 132 



^— e 



(£-*)■ 



2(c-2) 



£l-z>o(l^) 
z - > o (l^ 



In the intermediate regime where \l — £L\ -C L?- 1 , p{l) 
has the form 



p(Z) ~ L 



-c/(c-l) 



Vr. 



i - ^ 



LV(c-i) 



(61) 



Therefore p(Z) has a peak centered around Z ~ £L with a 
power-law decay on the right and a stretched exponential 
decay on the left. Integrating p(l) as given by (|6ip around 
Z « £L yields an order 1/L contribution implying the 
existence of a macroscopic loop. 

(b) For c > 3 the resulting behavior is summarized in 
Eqs.(lOO-lOl) of [1| and read 



-C) 2 

4b = < 0(L~ 1/3 ) 



Note that 62 < 0. Hence p(l) is of the form 



P(l) 



I ~ O(L) 
l\ « O (L 2 / 3 ) 



Therefore in this case the condensate bump has a Gaus- 
sian form with weight 0(1/L), as in the case 2 < c < 3. 

We conclude that the loop size distribution p(l) is a 
power law (reminiscent of the critical phase) for O(l) 
loops, superposed with a bump centered around £L = 
as shown in Figj3] The precise form of this con- 
densate peak differs from the one found in the grand- 
canonical analysis, although both ensembles yield the 
same phase diagram in the large L limit. This result is 
very similar to what is found in the context of ZRP, how- 
ever it is not exactly the same. Within the ZRP, above 
the critical density any further increase in the density is 
absorbed by the condensate. Here, on the other hand, 
the total length of the loops in the critical phase changes 
with temperature above T c . In particular, it is finite at 
T c and approaches zero at T — ¥ 00. Since the loop size 
distribution in the critical phase is fixed above the crit- 
ical temperature, it implies that the number of loops in 
the critical phase varies with T. 



V. CONCLUSIONS 



We analyzed the denaturation transition of circular 
DNA chains, assuming that opening denatured loops in- 
duces formation of supercoils. As in the case of non- 
circular DNA the thermodynamic behavior of the model 
is found to be determined by the loop entropy parameter 
c. We find that for c < 2 the model exhibits no transi- 
tion while for c > 2 the transition is continuous, of order 

■ Thus for c > 3 the transition is second order, 

while for 2 < c < 3 (which includes the physical value of 
c ~ 2.12) it is of higher order reaching oo-order as c — ¥ 2. 

In addition, the nature of the denaturated phase is 
rather different from that of the non-circular DNA. Here 
a macroscopic loop (condensate) is formed above T c 
whose length increases continuously as the temperature 
is increased. This is different from the denaturated phase 
in the non-circular case, where the two strands are fully 
separated at all temperatures above T c . This is remi- 
niscent of Bose-Einstein condensation and to similar real 
space condensation encountered in models such as the 
ZRP [33, HH. Furthermore, the difference observed in 
the condensate peaks of canonical and grand-canonical 
ensembles (for finite L) has the same mathematical struc- 
ture as in the ZRP. 

A different mechanism for absorbing the extra linking 
number produced by opening of loops in circular DNA 
has been considered previously [22I [23j . In this mecha- 
nism the extra linking number is compensated by over- 
twist of remaining bound segments of the molecule at 
the cost of an elastic energy. This mechanism also yields 
smoothening of the denaturation transition as obtained 
in the present paper. It would be of interest to consider 
the denaturation transition in the case where both over- 
twist and supercoils are present. 

Finally, our results apply to a homogeneous polymer 
where there is a single binding energy. It is well known 
that introducing disorder also smoothens the first-order 
transition in the PS model [33j. The influence of sequence 
inhomogeneity on the present melting transition which is 
already smoothencd by topological constraints is an open 
question. 

We thank O. Cohen, M.R. Evans, O. Hirschberg, S.N. 
Majumdar and E. Orlandini for helpful discussions. This 
work was supported by the Israel Science Foundation 
(ISF) and the Turkish Technological and Scientific Re- 
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